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ABSTRACT 

We examine the extent to which the self-annihilation of supcrsymmctric ncutralino 
dark matter, as well as light dark matter, influences the rate of heating, ionisation and 
Lyman-a pumping of interstellar hydrogen and helium and the extent to which this is 
manifested in the 21cm background signal. Unlike previous studies we fully consider 
the enhancements to the annihilation rate from dark matter halos and substructures 
within them. We find that the influence of such structures results in significant changes 
in the brightness temperature. The effect on the global signature at redshifts within 
the range probed LOFAR (i.e. z < 12) is on the edge of its sensitivity in the case 
of neutralino dark matter, and very likely to be detected for annihilating light dark 
matter. 
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1 INTRODUCTION 

The standard cosmological model, motivated by measure- 
ments of temperature anisotr opies in the Cosm i c Mi- 
crowave Background fCMB) (ppergel et all 120031 . 120071 ; 
iKomatsu" et al] 120081 ; iDunklev et al.l [20081) . the large scale 
distri bution of galaxies ( Cole et alj 20051 ; iTegmark et al.l 
120061 ) and by evidence of the accelerate d expansion of the 
Universe from Supernovae observations (|Astier et all [2006; 
IWood-Vasev et al 120071 ') . requires that the Universe is spa- 
tially flat, with a corresponding critical density, 30 percent 
of which consists of physical matter. However these obser- 
vations also indicate that only 4 percent of this matter is 
baryonic in nature, implying that the remaining 26 percent 
consists of an elusive, non-baryonic component called dark 
matter, because of the severe constraints that current astro- 
nomical data set on its radiative capabilities. 

Despite this compelling evidence for the existence of 
dark matter, its precise nature is still a topic of fierce de- 
bate. Particle physicists have independently supported dark 
matter by postulating the existence of a variety of exotic par- 
ticles with wide-ranging properties which may potentially 
solve problems in particle physics whilst resulting in a relic 
particle density consistent with current observational con- 
straints. 

The most intensely studied dark matter candidate is the 
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lightest neutralino (|Bertone. Hooper fc S ilk 2005), a weakly- 
interacting massive particle (WIMP) motivated by super- 
symmetric extensions of the Standard Model of particle 
physcs. In many of these extensions the neutralino is the 
lightest supersymmetric particle (LSP). In theories where 
the LSP is stable, for exam ple theories where R-Parity is a 
conserved quantum number (IWeinberJl 982; Hal l fc Suzukitl 



1 1984 lAllanach. Pedes fc Dreinerl Il999f ).~ the neutralino is 
thus a highly-motivated dark matter candidate. Further, an 
attractive feature of neutralinos is that a large region of 
the relevant supersymmetric parameter space can be inves- 
tigated using CERN's Large Hadron Collider (LHC), sched- 
uled for activation in 200qj. 

Whilst neutralino dark matter is 'cold', owing to its 
negligible free-streaming length (i.e. the length scale below 
which fluctuations in dark matter density are suppressed), 
warm dark matter (WDM) is typically lighter and possesses 
a much longer free-streaming length. WDM is a viable 
alternative to cold dark matter (CDM) models which may 
potentially resolve several shortfalls of the standard CDM 
model, like for example the over-predict ion of low mass satel- 
lites a nd the existence of cuspy halos (iHogan fc Dalcantonl 
20001; iDalcanton fc Hoganl l200ll; lAvila-Reese et alj I200T 



Colin. Valenzuela. fc Avila-Reesej 120081). Among WDM 



cand idates t here are sterile neutrinos (iDodelson fc Widrowi 
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iBerezinskv fc Valid Il993l ; lLattanzi fc Valid 1200?]) an d 
light dark matter (LDM) particles (|Boehm fc Favetll2004l) . 
What makes LDM interesting for this study is the fact 
that it can self-annihilate, as opposed to other forms of 
WDM, and therefore its annihilation rate can be enhanced 
by overdensities. 

In this paper we re-examine the influence of neutralino 
and LDM annihilations on the thermal history of the Uni- 
verse at times between the epochs of recombination and 
reionisation, referred to as the "Dark Ages", when gas ex- 
isted in a nearly uniform, dark, neutral state. The investi- 
gation of the Dark Ages is one of the frontiers of modern 
cosmology, and will be carried on by a new generation of ra- 
dio interferometers, such as LOFAE0, MWA0, 21 CMA0, and 
SKA0 These experiments will look for the redshifted 21 cm 
signal associated with the hyperfine triplet-singlet transition 
of neutral hydrogen. If dark matter annihilates or decays, 
the resulting products subsequently collide and heat sur- 
rounding gas, increasing its kinetic temperature and ionised 
fraction, which in turn is manifested as distinct features in 
the 21 cm background signal th at can be used to constrain 
the properties of dark matter dFurlanetto. Oh fc Pierpaolil 



trie properties ot dark matter (furlanetto, Un & fic 
12009 : IShchekinov fc Vasilievll2007l : IValdes et alJl2007l) 

W ith the exception of the very recent work bv lChuzhovl 
|2007l ). past studies proclaim that the heating effects as- 
sociated with the annihilation of WIMP dark matter are 
too small to be detected by current radio interferometers. 
However, these studies overlook the enhancements to the 
dark matter a nnihilation rate in ga lactic halos and in their 
substructures l|Tavlor fc Silkl l2003h . which could be large 
enough to make the dark matter signature detectable by 
the next generation radio telescopes. In this paper we calcu- 
late the effect of neutralino and LDM annihilations on the 
21 cm signal when accounting for the effect of dark matter 
clustering. 

The rest of this paper is organized as follows. In §[2] we 
elaborate on the basic properties of neutralinos and LDM. 
We also discuss the basic physics describing the way in which 
energy from annihilations is injected into the intergalactic 
medium (IGM) . In §0 and §2] we calculate the enhance- 
ment in the dark matter annihilation rate caused by the 
presence of halos and their substructures. In §[5]we estimate 
how much of the energy produced in a single dark matter an- 
nihilation is actually injected into the IGM. In §|B]we discuss 
the modifications to the differential equations describing the 
evolution of the ionised fraction and kinetic temperature of 
the IGM and subsequently use these equations to calculate 
the modified 21cm background. In §[7] we calculate the pre- 
dicted 21 cm background for our benchmark neutralino and 
LDM models and compare the results to the expected sen- 
sitivity of LOFAR, in order to assess the potential for a 
detection. Finally, in §[5] we summarise our results and draw 
our conclusions. 



2 DARK MATTER CANDIDATES 

The lightest supersymmetric (SUSY) neutralino is a super- 
position of higgsinos, winos and binos. Consequently neu- 
tralinos are electrically neutral and colourless, only inter- 
acting weakly and gravitationally and making them very 
difficult to detect directly. In SUSY mode l s that conserve 
R-pa x ity, the LSP is stable (|Weinberdll982l ; lHall fc Suzukitl 
1 1984 lAllanach et~aH Il999h . Consequently, in a scenario 
where present-day CDM exists as a result of thermal-freeze 
out, the dominant species of CDM could quite possibly in- 
clude the LSP. The relic density of the LSP will then heavily 
depend on its mass and annihilation cross-section. Through- 
out this paper we assume that the LSP is the lightest SUSY 
neutralino. The neutralino is a popular candidate for CDM 
because the theoretically-motivated values of these parame- 
ters correspond to relic densities which are of the same order 
as current measurements for f2cDM (for a more detailed re- 
view of the various p roperties and motivat ions for neutralino 
dark matter see e.g. iBertone et ail (|2005l )). 

Neutralinos possess a wide-range of annihilation spec- 
tra owing to the vast extent of currently available SUSY pa- 
rameter space. If the neutralino is lighter than the W and 
Z bosons, annihilations will be dominated by the process 
XX ~ * bb with a minor contribution by XX ~ * t + t~ . Assum- 
ing annihilations are dominated by the former process, the 
resulting spectrum will depend entirely on the LSP mass. 
For heavier LSPs, the annihilation products become more 
complex, often determined by several dominant annihilation 
modes including XX ~~ > W + W~ , XX ~ > ZZ or XX ~ * & as 
well as XX ~ * bb and XX ~ * t + t~ . 

The other dark matter candidate we are considering 
is LDM, consisting of MeV mass particles, which anni- 
hilate to electron-positron pairs and consequently were 
considered to be a possible source of the positrons con- 
tributing to the 511 keV positronium decay signature from 
the centre of the gala xy observed by SPI /INTEGRAL 
(Knodlsedcr et al. 2 0051 ) . Relevant analyses of this emission 
impose the constraint on the LDM mass kidm < 20 MeV in 
order not to overproduce de tectable gamma-rays from inner 
bremsstrahlun g processes (iBeacom, Bell fc Bertond 120051 ) 
(although see iBoehm fc Uwerf l 20061 )). A stronger, albeit 
less conservative constraint, tudm < 3 MeV can be ob- 
tained if one considers the generation of gamma-rays from 
the in-flight annihilation between positrons produced from 
LDM annihilation and electrons residing in th e interstellar 
medium of our galaxy (|Beacom fc Yuskel|[2003l ). 

Both in the case of neutralinos and LDM, the average 
rate of energy absorption per hydrogen atom in the IGM at 
a redshift z is given by 



1 Tl 

£ = n/abs PM '° (cTannf )m dm (l + zf C(l + z) (1) 

2 n H ,o 

where rriDM is the mass of the dark matter particle, (iTann.w) 
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6 MeV LDM particles can also potentially annihilate directly 
into neutrinos and photons. However most theories suppress this 
emission in order to be consistent with observational constraints. 
Here we only consider scenarios where LDM annihilates entirely 
to electron positron pairs, so that our results can be considered 
as an upper limit to the more general case. 
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is the thermally-averaged dark matter annihilation cross- 
section, 71dm, o and nn,o are the current average number den- 
sities of dark matter and hydrogen respectively, and / a bs. 
is the fraction of energy which is absorbed by the IGM. 
C(l + z) is the redshift-dependent enhancement of the an- 
nihilation rate owing to dark matter structures, relative to 
a completely homogeneous Universe. Q 



3 EXTRAGALACTIC DARK MATTER 
ANNIHILATION RATE 

In the standard cosmological model all structure in the uni- 
verse originated from small amplitude quantum fluctuations 
during an epoch of inflationary expansion shortly after the 
Big Bang. The linear growth of the resulting density fluctu- 
ations is then completely determined by their initial power 
spectrum, which for ACDM is usually assumed to be a power 
law with spectral index n. Current limits on n from ob- 
servations of temperature fluctuations in the Cosmic Mi- 
crowave Background conducted by the Wilki nson Microwave 
Anisotropy Probe, hwmap = 0.963lo;oi5 (|Dunklev et al.l 
120081 ). supports the existence of a power spectrum consis- 
tent with inflation. 

During the expansion of the Universe, the small ini- 
tial density fluctuations will eventually grow and produce 
the structures we observe today. In the currently accepted 
cosmological model, smaller structures form first and then 
merge to form larger ones in a process of "bottom-up" hier- 
archical structure formation. The mass distribution at any 
given redshift can potentially be determined through the use 
of numerical simulations. 

As a first approximation, the smaller progenitors form- 
ing larger isolated structures are completely disrupted after 
merging and the resulting 'smooth' dark matter density dis- 
tribution can be described by a continuous function conven- 
tionally of the form 



p(r) 



(r/r s y< [l + (r/r fl )"] ( ' 3 - 7)/o 



(2) 



where r is the distance from the centre of the halo (spheri- 
cally symmetrical) , r s is a scale radius, p 3 is a normalisation 
factor, and a, j3 and 7 are free parameters. 

However, N-body simulations of CDM halos reveal that 
a wealth of substructure hal os (subhalos hencefort h) exists 
within such halos. Moreover, iDiemand et al.l (|2008h recently 
claimed, based on the results of the second generation of 
Via Lactea simulations, that a further generation of sub- 
subhalos exists with a self-similar mass distribution relative 
to the parent subhalo. This suggests the possibility that if 
one were to conduct simulations with sufficiently high res- 
olution, one would find a long nested self-similar series of 
halos within halos within halos etc., all the way down to the 
smallest halos. 

This intriguing result has significant implications for the 
indirect detection of annihilating dark matter, since the rate 
of dark matter annihilations is proportional to the square of 
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the local density, hence the presence of over-densities can 
significantly increase the annihilation rate relative to that 
obtained with a smooth background distribution. 

The above scenario applies to structures formed in a 
CDM-dominated Universe. In a WDM-dominated Universe, 
the significant damping of small-scale density fluctuations, 
due to the larger free-s treaming length, sho uld be taken 
into account. Following iBardeen et al.l (|l986t ) this can be 
accounted for by using the modified power spectrum 

P{k) = T^ dm (/c)P A cdm, (3) 
where the WDM transfer function is approximated by 



Twdm(A:) = exp 



kR f _ (kR f ) 2 
2 2 



(4) 



where Rf is the free-streaming length. 

For WDM particles with negligible interaction rates the 
free-streaming length is related to the particle mass didm by 
jBardeen et al.lll986l ) 



R f , n = 7.4 x 10" 



h 



lMeV 



0.258 



1/3 



0.719 



5/3 



hT 1 Mpc. 



(5) 



However, the interaction rates for self-annihilating LDM in 
the models we consider are non- negligible. In this case, the 
free-streaming length is given by (|Boehm fc SchaefferlliooBT l 



Rf.i 



0.3 



6X lO^S-Hl + Zdec) 3 

lMeV\ 1/2 „ 
Mpc, 

TTlDM / 



1/2 



(6) 



where Fd cc ,DM is the WDM self-annihilation rate at the de- 
coupling redshift Zdec. 



1 Pc,0^DM,0 



TTlDM 



(Cann.U )dcc.(l + ^dec) 



(7) 



where (<T a nn.w)dec. i s th e thermally-averaged product of the 
WDM annihilation cross-section multiplied by relative speed 
evaluated at the Zd cc , which in order to obtain the cor- 
rect relic density today is required to be approximately 
(oWi.u}dec. ~ 10~ 26 cm 3 s _1 . 

For 771dm = 3MeV, R f . n = 2.38 pc and R f<i = 97.99 pc, 
while for ttidm = 20MeV, Rf >n = 0.190 pc and Rf ti = 
14.70 pc. Hence, in both cases the comoving free-streaming 
length set by interactions is at least an order of magnitude 
larger than that when interactions are completely negligible, 
and consequently we must invoke the former in our deter- 
mination of the cut-off scale in the WDM power spectrum. 

We follow lAvila-Reese et all (120011 ) and define a 
characteristic free-streaming wavenumber kf such that 
TwDM(fc/) — 0.5 leading to kf ~ 0.46/7?/. This wavenumber 
is then related to a characteristic filtering mass Mf by 



47T _ f Xf 

M f = — pwdm I — 



(8) 
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where Xf = 2n/kf = 13.6-R/. In this paper we invoke the ap- 
proximation Af m in. ~ Mf, where here M m i n . is the minimum 
mass of a LDM halo, and equal to approximately 46 Mq and 
0.16 M Q for m dm = 3MeV and 20MeV respectively. Since 
the mass within a given comoving volume is constant as the 
Universe expands, the result is independent of redshift. 

Below we perform a series of detailed calculations illus- 
trating the enhancement of the annihilation rate relative to 
that obtained with a completely smooth Universe, known as 
the clumping factor. 



4 CALCULATION OF THE CLUMPING 
FACTOR 

We assume a standard homogeneous, isotropic Universe with 
a flat spatial geometry. Let R(M, z) be the average annihi- 
lation rate per baryon in a generic dark matter halo of mass 
M located at redshift z. Even for large M, this source can 
be regarded as an unresolved point-source and we assume 
this throughout for all halos considered. The rate of annihi- 
lations at a given redshift is then equal to 



T(z) = (l + z) 



dM—(M,z)R(M,z) (9) 



where we have introduced the unconditional halo mass func- 
tion, dn/dM, i.e. the comoving number density of virialised 
halos with mass M located at redshift z, (the factor (l + z)' 5 
converts this from comoving to proper density). The inte- 
gral spans over the mass range M > M m - m ., where M m i n . 
can be as small a s 10~ 12 Mq due to k inetic decoupling in 
the case of CDM (jProfumo et al.l[2006T ), and approximated 
by the filtering mass (JH| in the case of WDM. 

Three ingredients are required in order to calculate the 
annihilation rate ([9}. We need to specify the annihilation 
cross-section of our dark matter candidates (in our case neu- 
tralinos or LDM). Secondly, we need to specify the dark mat- 
ter density profile of a generic halo of mass M at redshift z. 
Finally we need an estimate of the distribution of halos, i.e. 
an estimate of the halo mass function dn(M, z)/dM. 



4.1 The halo mass function 

Press-Schechter theory (|Press fc Schechter]|l974l ) postulates 
that the cosmological mass function of dark matter halos 
can be expressed in the universal form 



dn _ _po_ ,, s dlogz/ 
dM ~~ M 2 " 1 ^' dlogM 



(10) 



where po is the average comoving dark matter density pa = 
PcQm, where p c is the present critical density of the Uni- 
verse. The parameter v — S sc /a(M), is defined as the ratio 
of the critical overdensity required for spherical collapse at 
redshift z, extrapolated using linear theory to present time, 
and a(M), the r.m.s. of primordial density fluctuations when 
smoothed on a scale which contains mass M, again extrapo- 
lated using linea r theory to present tim e. The form of S sc (z) 
can be found in iTegmark et ah! (|2005l ). a(M) is related to 



the power spectrum P(k) of the linear density field extrap- 
olated to present time by 



a 2 (M) 



J d :i kW 2 (kR)P(k), 



(11) 



where W is the top-hat window function at the scale R = 
(3M/47rp) 1//3 where p is the mean matter de nsity. We utilise 
the a nalytical approximation specified in ITegmark et al.l 
I2006T ). relevant in the linear regime long after the relevant 
fluctuation modes have entered the horizon, when all modes 
grow at the same rate, which means that a(M) can be fac- 
tored as a product of two functions, one solely dependent on 
redshift z and the other solely dependent on the comoving 
spatial scale R. We normalise P and a by computing a at 
R — 8ft _1 Mpc and setting the result equal to the cosmologi- 
cal parameter as as m easured by WMAP, as = 0.796±0.036 
(|Dunklev et aLlbOOST ). 

The fir st-crossing distribution f(v) has the following an- 
alytical fit (|Sheth fc Tormenl2002T ) to the N-body simula tion 
results from the Virgo consortium (|jenkins et al.|[T998l ) 

vf(v) = A[l + (av)-*>] (g) 1/2 ex p(-^) (12) 

where the parameters p = 0.3, and A = 0.32218 are deter- 
mined by the requirement that all mass lies within a given 
halo, i.e. J dvf(v) — 1 or equivalently J dMMdn/dM — po. 



4.2 The density profile of dark matter halos 

Since the rate of dark matter annihilation scales with den- 
sity squared, it sensitively depends on the density profile 
of each halo. We consider three universal density profiles 
to model the smooth distribution of dark matter within 
each halo (substructure will be dealt with later in this sec- 
tion). Firstly we conside r the popul ar profile proposed by 
iNavarro. Frenk fc White! (|l996l . ll997T ) (NFW), 



p(r) = 



(13) 



(r/r s ) [1 + (r/r.)]' 

which corresponds to a = 1, j3 = 3 and 7 = 1 in eq.Q. 
Secondly, we consider a profile with a significantly large r 
slope, specifically the one proposed by iMoore et al.1 (| 19991 ). 



p(r) 



Ps 



(r/r,) 1 - 5 [1 + (r/r a ys\ 



(14) 



which corresponds to a = 1.5, (3 — 3 and 7 = 1.5. Both of 
these profiles have the same functional form and are both 
singular towards the galactic centre (in fact, the slope of the 
Moore profile must necessarily be truncated for r < r m i n ., 
where r m i n . ~ (see below), otherwise the integral of density 
squared diverges. 

However, there have been indications that cuspy profiles 
are inconsistent with observations, specifi cally regarding the 
rot at io n curve s of s m all-scale galaxies dFlores fc Primackl 
1994 IMoore] Il994l; [Weldrake de Blok fc Walter! 120031 ; 



Donato. Gentile fc Saluccill2004l ; iGentile et al.ll2007r i. which 
are more likely to be consistent with density profiles possess- 
ing flattened cores similar to that which may be achieve d 
with WDM (|Hogan fc Dalcantonl I2OO0I ; IColin et al.l 120081 ). 
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Therefore, we lastly consider the Burkert density profile 
(|Burkertll 19951 ) 



p{r) = 



[1 + (r/r 3 )] [1 + (r/r s )2] 



(15) 



which has been shown to be fairly consistent with 
the rotation curves of a large number spiral galaxies 
(ISalucci fc Burkert|[2000h . 



4.4 Clumping factor for smooth halos 

We are now in a position to be able to calculate the clumping 
factor C(l + z) attributed to extragalactic halos with smooth 
density profiles described by equations (|13[) -(115 [ ), with con- 
centrations given by equation (|19l l. 

We start by writing down the annihilation rate R(M, z) 
within a dark matter halo of mass M located at redshift z, 
given by 



4.3 Concentration - mass relation for dark matter 
halos 

Here we introduce the virial concentration parameter c v i r ., 
defined by 



^*vir. 



(16) 



where r 3 is the scale radius defined above and r v i r . is the 
virial radius of the halo, defined to be the radius encapsu- 
lating the virial mass M of the halo within which the average 
density is equal to the overdensity A v i r . times the average 
cosmological density p(z) at that redshift 



(17) 



For A yjr. we use the approximation given in lTegmark et all 
(1200611 



A vir . ~ 18tt 2 + 52.8a; - 7 + 16a;, 



(18) 



where x(z) — £l\(z) /Q.m(z), (A v i r . ~ 311 at z = for Qm = 
0.3 ad Qa = 0.7) and accurate to within 4% of the exact 
numerical calculation at relevant times. 

There has been evidence from simulations revealing a 
strong correlation between the halo mass M and its cor- 
responding concentration c v i r ., with larger concentrations 
in smaller mass halos, which is consistent with the idea of 
bottom-up hierarchical structure formation with smaller ha- 
los collapsing at earlier time s when the average density of the 
Universe was much greater (iNavarro et al.lll996l.ll997l ) . This 
relationship was later reaffirmed by Bullock et al.l (|200ll ) 
(B2001 hereafter) using a sample of simulated halos in the 
mass range 10 11 < M/h -1 Mq < 10 14 , who proposed a toy 
model to describe this behaviour, which is popular in the rel- 
evant literature: On average, a collapse redshift z c is assigned 
to each halo of mass M through the relation M, = FM, 
where at a redshift z the typical collapsing mass M*(z) is 
defined implicitly by the relation <j(M*(z)) = S sc (z) and 
is postulated to be a fixed fraction F of M, which follow- 
ing ImhoeFai] (|2002[ ) we set equal to 0.015. The density of 
the Universe at redshift z c is then associated with a char- 
acteristic density of the halo at z. Therefore, here we use 
the average concentration-mass relation obtained using the 
above method, which is given by 



.(M,z) = K 



l+_gc 
1+Z 



where K ~ 5, for f2 A = 0.742, Q , M 
a 8 = 0.796 jDunklev et al.ll200Sl ). 



Cvir.(M, z = 0) 

1 + z 
0.258, h = 



(19) 
0.719 and 



R(M,z) = -— 2 — 



2 mJ M Tib(z) 



p 2 (r)4nr 2 dr. (20) 



where {a&nn.v) is the thermally-averaged product of the 
dark matter annihilation cross-section multiplied by relative 
speed, and n^z) is the average cosmological baryon density 
at redshift z. For a halo of mass M with an NFW profile 
with concentration c v j r .(M, z), R(M,z) is easily calculated 
to be 



R{M,z) 



1 ((Tann.w) / p s \ 4lT fr vil .(z,M) 



2 n b (z) \m D M J 3 \c vir . (z, M) 
1 



= xU- 



[l + c vit .(M,z)Y 



(21) 



By equating (|17p for the virial mass M to the integral 



A I 



J p(r,c vir .(M, z))47rr 2 dr 



4tt 



Ps(M,z) 



c vir . (M, z) 
log(l + c vil Xz,M)) 



(l + c vir .(z,M)) 



(22) 



we obtain the relation for the scale density 
M 



p s (M,z) 



4tt 



(c vir V (M, Z )) 



1 



[iog(i+ Cvir .(^M)) - 

(23) 

For the Burkert profile, it is not possible to express the re- 
sults corresponding to equations (1 2 2 ll and (|23p in analytical 
form. For the Moore profile, in order for the integral over 
density squared to be finite we must truncate the density 
below a radius r m i n .. Defining the variable x = rc v i r ./V v ir. 
with Train. = rmm.Cvi r ./ r vir., expressions corresponding to 
((22)) and ((23j) for the Moore profile 



for r > r min . 



P(r) (r/r 3 )i' 5 [1 + (r/r,) 1 ' 5 ] 

p(r) = p(r min .) for r < r min 

are calculated as 



(24) 
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R(M,z) = 



1 <(7a 



2 4tt /r vir .(,z,M) N 3 



2 ra 6 (.z) \m D in/ 3 \Cvir.(z,M) 

X F\ (Cvir. i ^min.)] 



(25) 



where 



F\ (Cvir. ) 3*min 



+ 



loji 



3(l + z min .) 2 
1 

L5 
1 

L5 



^L(l+4ir.) 



1 + ^ 



(26) 



and 



p s (M,2) = 



M 



(c vir V (M, 2 )) 



4tt 



1.5 



rloe 



1 + 

M 



1 a; 1 ? 



3(l + ^f n .) 



(c vir V (M,2)) -^ 2 ( Cvir - ' x min.) 



4tt 



(27) 



respectively. 

Then following equation © for the contribution to the 
dark matter annihilation rate per baryon Fhaios(z) by halos 
located at redshift z, we obtain 



(W maXi 

j dM^(M,z)R(M,z) 



1 (O'ann.l') 

2 m^ M n b (z) 



(1 + *) 



■.<M,z) 



(28) 

To calculate r m i n . for the Moore profile we adopt the 
following method: Within some minimum radius we assume 
that most of the neutralino dark matter that we are consid- 
ering would have self-annihilated, leaving a flattened density 
core. The size of the core will depend on several time-scales: 
the time-scale for a single dark matter annihilation (per par- 
ticle) tann. ~ (n x (av})~ , the average time-scale for the re- 
plenishment of an annihilated dark matter particle within 
the core owing to in- falling dark matter from larger radii tin. , 
and the collapse time-scale of the halo in question t co n <C th, 
where th ~ 10 17 s is the Hubble time. Clearly, for a core to be 
produced the average time-scale for replenishment of a dark 
matter particle within the core must be more than the time 
it takes to annihilate, and must be much smaller than the 
time-scale for the formation of the halo, which if the halo is 
to exist must be much smaller than the Hubble time. Hence 
for a core to form we require the condition 



. , 1 1 1 

n«H > 7— » 7 — » 7- 

Wn. ^ coll. 



(29) 



Here, we adopt the conservative criterion 
(Px( r min.)/m x )(cru) ~ tj~ , with canonical values of 
the neutralino mass and annihilation cross-section of 
m x ~ 100 GeV and (av) ~ 10~ 26 cm 3 s~\ Since the Moore 
profile is monotonically decreasing with increasing radius, 
this sets an upper limit for x m - m . which is approximately 
equal to 10~ 8 for the Milky Way at present day, which is 
consiste nt with simila r app roximations by other authors 
(e.g. see lTavlor fc Silkl (120031 )1. 

The corresponding rate of dark matter annihilation per 
unit volume contributed by the smooth background density 
at redshift z is given by 



r , , _ 1 (o-ann.u) _2 , v 

J- smooth^; — r — _ , \ PDM( z ) 

2 rib(z) 



(30) 



where Pdm(z) = p c ,o^DM,o(l + z) 3 . Therefore, we define the 
clumping factor for smooth halos Chaio(z) as 



Chalo(z) = 1 + 
= 1 + 



r h alo(^) 



Tsmooth (■£■) 

(i + ^) 3 



Pdm( 2 ) 

••vir-CM,.) 

J p 2 (r)47i 



(31) 



where Ch a io(«) — > 1 for a completely smooth Universe. 
For halos with NFW profiles, using equations (|21fl . (|23fl . 
and (1301) we obtain 



Chalo(z) = 1 + 
(1 + Z) 3 

pIm( z ) 



z)p 2 a (M,z) 



4tt / r vll .(M,z) 



3 \c vil .(M,z) 



(l + c vir .{M,z)Y 



(32) 



whilst for halos possessing a Moore profile, using equations 
((25)) , l|26)l. (f27|). (p8)) and (O we obtain 



Cnalo(«) = 1 + 

(1 + ^) 3 
Pdm( z ) 



/ 



kf-^(M,z)p 2 (M,;z) 



dM 



4tt /r vir .(M,z) 



3 V c vir.(M,z) 



F , i(c vir .(M,z),2: min .(M,z)). 



(33) 



Since the integrals over powers of the Burkert density profile 
are non-analytical we do not attempt to simplify its corre- 
sponding expression for Chalo beyond equation (pi) . 

In figures [TJ3] we display plots of Ch a io(z) as a function 
of z for halos with (a) NFW profiles, (b) Moore profiles and 
(c) Burkert profiles. 
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Figure 1. Clumping factor as a function of rcdshift for dark 
matter halos with mass M > Af m j n . with smooth NFW density 
profiles with a c v ; r — M relation described by eq. 1191 1, trun- 
cated at a halo mass M cu t. The displayed curves correspond to 
values of (A/ min ./M Q ,M cu t/M ) equal to (10~ 12 , 10~ 12 ) (thin 
black solid), (10- 12 ,10 6 ) (thin blue dashed), (10- 4 ,10~ 4 ) (thin 
red dot-dashed), (10~ 4 , 10 6 ) (thick green dashed),and (10 6 ,10 6 ) 
(thick magenta dot-dashed). 



Halos with cuspy density profiles, such as the NFW and 
Moore profiles, are typical of CDM halos for which the min- 
imum mass cut-off scale in the matter power spectrum is 
determined by collisional damping and free streaming in the 
early Universe. For WIMP dark matter the value of M m i n . 
can range from 1O~ 12 M0 to 1O _4 M0 for typical kinetic 
decoupling temperatures. Hence in Figs. [1] and [2] we illus- 
trate the effect on the clumping factor for values of M m i n . of 
10~ 12 , 10~ 4 and 1O 6 M0, where the latter value is the typical 
minimum mass of resolved subhalos in simulations of galac- 
tic halos. We also demonstrate the influence of truncating 
the Bullock et al. concentration-mass relation (referred to 
as the "B2001 relation" hereafter) below a mass of 1O 6 M , 
aswell as using the relation when extrapolated to M m i n .. 

In Fig. [3] we show the clumping factor for halos ha- 
los with flattened cores like the ones possibly formed by 
WDM. In particular, we plot Chaio for minimum halo masses 
M m i n . — 46 M© and O.I6M0, corresponding to the values of 
the damping mass |(HJ obtained using ttiwdm = 3MeV and 
mwDM = 20MeV, respectively. We also again illustrate the 
effect of using relation eq. (|19|l when extrapolated to M m i n . 
or truncated at 1O 6 M0. The selected values of mwDM cor- 
respond to the respective upper limits on the LDM particle 
mass from constraints on in ner bremsstrahlung g amma ray 
flu x from the galact i c cent re (|Beacom et al ] |2005t ) (although 
see iBoehm fc Uwerl (I2006T I). and from in-flight annihilation 
(|Beacom fc Yuskellbooi T between positrons from LDM an- 
nihilation and electrons in the interstellar medium. 



4.5 Clumping factor for halos possessing 
sub-halos and sub-sub halos 

Thus far we have considered the amplification of the dark 
matter annihilation rate for isolated halos with smooth 
density profiles. However, as already mentioned, N-body 
simulations indicate that a significant proportion of the 
smaller progenitors giving rise to larger mass halos survive 
the merging processes and the tidal forces exerted upon 
them during their orbital motion within halos. In partic- 
ular, the most recent ACDM Via Lactea II simula tions 
of galactic halos presented in iDiemand et al.l (|2008l ) and 
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Figure 2. Same as for figure[T] except now for halos possessing 
Moore density profiles. 




Figure 3. Same as for figure[TJ except now for halos possess- 
ing Burkert density profiles for values of (M m ; n /Mq, M c ^/Mq) 
equal to (0.16,0.16) (thin black solid), (0.16, 10 6 ) (thin blue 
dashed), (46,46) (thin red dot-dashed), (46, 10 6 ) (thick green 
dashed) and (10 6 , 10 6 ) (thick magenta dot-dashed). 



in lKuhlen, Diemand fc Madaul (|2008l ) (KDM hereafter), re- 
vealed a second generation of surviving substructures within 
halos (designated as "sub-subhalos" ) . Further, these simu- 
lations suggest that the mass distribution of sub-subhalos 
within their host subhalo is the same as the mass distribu- 
tion of subhalos within their host halo. 

Since the dark matter annihilation rate scales with den- 
sity squared, such subhalos and sub-subhalos could provide 
significant enhancement to the annihilation rate even for 
modest substructure mass fractions within halos/subhalo. 
For halos of mass M these have been suggested to be as 
much as 10% fo r subhalo masses M s in the range 1(T 5 < 
M s /M < 1CT 2 (|Diemand et al.1 |200S( ) (which corresponds 
approximately to a constant mass fraction per subhalo mass 
decade of 3%, owing to the fact that the subhalo mass func- 
tion has a slope of approximately 2, see below). However, 
owing to the fact that substructures invariably form earlier 
than their host halos, and that tidal disruption is unlikely to 
effect the inner density profiles of structures (i.e. where the 
majority of the enhancement originates), the concentration 
of substructures may be significantly greater than that of 
their host halos. The simulation results recently presented in 
KDM are consistent with the ratio iV c = c^°/c^ bhal0 ~ 3 
for subhalos located at solar radii within galactic halofjf], 



8 Although N c demonstrates a slight galactocentric radial depen- 
dence, the authors of KDM claim that the effect on the overall 
annihilation rate is negligible. 
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Figure 4. Plots of the total clumping factor from halos and sub- 
halos, given by Ctotal.O) = 1 + Chalos(z) + C su bhalos(«), where 
we have renormalised the contribution from halos to account for 
their finite mass fraction in the form of subhalos per subhalo mass 
decade F su b. for varying values of the substructure mass function 
index (3. Displayed are values of C tota i for structures with NFW 
(purple) and Moore (blue) density profiles for (3 = 2 and 1.8 (solid 
and dashed curves respectively) and (3 = 2 and 1.6 for the Burk- 
ert profile (black solid and dashed curves respectively). We take 
M cut = Af min . = lCr 12 Af Q for NFW and Moore profiles and 
II. Hi M . for Burkert profiles. In all cases -F su b. = 3%, N c = 3. 



where A is the appropriate normalisation of AN/AM S , p s s u ' 
is the subhalo scale density, which can be obtained from 
the previous expressions (|23[) or (|27p for NFW and Moore 



profiles respectively with the substitutions c v ir 



and 



M — > M a . Then integrating this contribution over all halos 
at redshift z we obtain the annihilation rate for all subhalos 
residing within such halos 



Tsubhalos (^) 



M MM : z) R suh .(M,z,F sub ) 



AM 



(o-ann.w) 



2m 2 

DM 



x / AM 



An(M,z] 



AM 

r vil .(z,M s ) 



A(M,F Eub .) J 



am s m: 



P {r,c v ^([M s , z])4-Kr Ar, 



r=0 



(36) 



and following equation (|31|) . we obtain the associated sub- 
halo clumping factor 



whilst the numerical simulations of Bullock el al. show that 
on average iV c ~ 1.5 for halos of mass M ~ 5 x 10 11 M Q 
(B2001). 

Here we calculate the contribution to the clumping fac- 
tor by halos possessing substructures with a self-similar mass 
distribution. Consider a dark matter halo of mass M with a 
subhalo mass distribution function given by 

AN(M) 



AM 3 



oc M, 



(34) 



where the index (3 is assumed to be time-independent and 
approximately equal to 2, i.e. equal mass per decade in sub- 
halos (KDM). 

There are indications that (3 may slightly deviate 
from this value, particularly for WDM substructures. 
iKnebe et al.l (|2008l ) claim that (3 may be as small as 1.6. 
However, as shown in Fig. [4] the effect on the clumping 
factor of varying (3 deviating slightly from 2 is small. We 
therefore adopt the value (3 = 2 for both CDM and WDM. 
Consequently, each subhalo mass decade contributes a con- 
stant fraction F Bub , of the halo mass. 

Adopting a course of reasoning analogous to that used 
to derive equation (|20|l . the rate of annihilation within a 
similar halo, owing solely to subhalos within it, possessing 
smooth density profiles p(r), is then given by 



R sub .{M,z) = 



((Taan.v) f ^ AN{M,F sub ,) 



2m 2 

r vir . (z,M s ) 



AM a 



p 2 (r,^t;[M s ,z])Anr 2 Ar 



r=0 

Im 2 

r vir (z,M s ) 



A(M,F suh .) / AM 3 M, 



J p 2 (r,^;[M s ,z])4* 



(35) 



Csubhalos — 1 H - 



Tsubhalos ( -2- ) 



Tsmooth (^) 



= 1 + %±f^ / dM^^l A(M,F sub .) f AM S M- 
PdmW J AM J 

r vir . (z,M a ) 

x j p 2 (r,clt-[M s ,z])47 



iirr 2 Ar. 



(37) 

However as mentioned above, each subhalo is likely to 
itself host substructures with mass function 



diV 
AM as 



A(M s ,F sa ,(3)M- 



(38) 



Owing to the self-similar nature of the mass distribution 
of substructures within halos, we take the values of the in- 
dex (3 SS and the sub-subhalo mass fraction per mass decade 
F as to be equal to (3 and F sub . respectively. Hence, follow- 
ing the above treatment for halos and their subhalos, the 
clumping factor for all sub-subhalos with virial concentra- 
tion Cvi r . residing within subhalos, themselves residing within 
halos located at redshift z is given by 



O su b-subhalos = 1 + — TT I - 



x A(M S ,F S ) f AM S M7 P 



x A(M 3S ,F S3 ) J AM SS M~ 



/T vlr . (z,Ms) 

/ p 2 (r,c v s ir .[M ss , z])4irr 2 Ar. 

Jr=0 



(39) 

Finally, using equations (|31|l , (|37p and (|39[) , we can now 
write down the total clumping factor for all structures at 
redshift z: 
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Figure 5. Plots of C to ta,l( z ) f° r structures with Burkcrt (solid 
black), NFW (dashed black) and Moore (dot-dashed black) den- 
sity profiles for values of (A/ min ./A/ Q , A/ Cut /A/ Q ) = (10 6 ,10 6 ) 
with subhalo and sub-subhalo mass fractions per decade F su b 
and F ss of 0.3 and a relative concentration ratio N c of 3.0. Also 
displayed are the corresponding values of Ch a io( 2 )i C au bhalos( 2 ) 
and C BU b_ su bhalos( 2; ) associated with those structures with NFW 
profiles (dashed blue, red and green curves, respectively). 
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Figure 6. Plots of C to tal( 2 ) for structures with Burkert (solid), 
NFW (dashed) and Moore (dot-dashed) density profiles for val- 
ues of (A/ min ./A/ , Afcut/M©) = (10 6 ,10 6 ) with subhalo and 
sub-subhalo mass fractions F su b and F BB of 10% and a relative 
concentration ratios N c of 3.0 (black curves) and 1.5 (magenta 
curves) . 



Ctotal = 1 + (ChaloO) - 1) 
+ (C su bhalos(z) — 1) 

+ (C S ub-subhalos(z) — 1), (40) 

where we must modify the normalisation of expressions Chalo 
and C su bhaios to include the fact that a proportion of the 
mass of each halo and subhalo is provided by substructures. 
From the repetitive structure of equation (f39|) . one can eas- 
ily see how to extend the present scenario to include higher 
generations of substructures, but considering there is no ev- 
idence for such structures we do not include them in this 
study. 

In figures (JSJ) , © and Q , we display the total clumping 
factor for various different scenarios in order to illustrate the 
relative ambiguity in its values associated with the varying 
levels of uncertainty on the various parameters involved. In 
addition in fig. ((5JI we illustrate the halo, subhalo and sub- 
subhalo contributions to the total clumping factor Ctotal (z) 
for structures possessing NFW density profiles. One can also 
see that the relative contribution to Ctotal (z) steadily de- 
scends for realistic values of the concentration ratio N c and 
substructure fraction F su b to the extent that further gener- 
ations of substructures, if they exist, are unlikely to increase 
Ctotal (z) by more than a few percent, and hence we neglect 
their contribution here. 



5 ABSORBED FRACTION OF INJECTED 
ENERGY 

In order to compute the rate at which energy is injected 
into the IGM, we need to reliably estimate how much of the 
energy produced in the DM annihilations is actually trans- 
ferred to the IGM, possibly altering its thermal/ionisation 
history. The energy absorbed fraction term / aDS in Eq. {TJ 
thus denotes the fraction of the particle rest mass energy 
that is injected in the IGM at a given redshift. This de- 
pends on the spectrum of DM annihilation products and on 
their interaction with the IGM. 




1+z 



Figure 7. Plots of Ctotal ( 2 ) for structures with Burkert (ma- 
genta solid), NFW (blue dashed) and Moore (black dot-dashed) 
density profiles with subhalo and sub-subhalo mass fractions 
per decade F su b and F ss of 0.03, a relative concentration ra- 
tio N c of 3.0, with values of (M min JMQ, M cu t/M Q ) equal to 
(0.16,0.16), (0.16, 10 6 ), (46, 46), (46, 10 6 ) and (10 6 ,10 6 ) for 
the Burkert profile (respective descending curves at high-2), 
and (10- 12 ,10- 12 ),(10- 12 ,10 6 ),(10- 4 ,10- 4 ),(10- 4 ,10 6 ) and 
(10 6 ,10 6 ) for NFW and Moore profiles (respective descending 
curves at high- z). 



Neutralino dark matter can annihilate directly into ei- 
ther a fermion pair or weak gauge bosons. Since the cross 
section for annihilation to fermion pairs is proportional 
to the square of the final state fermion mass, this pro- 
cess will be dominated by heavy final states, namely 66, 
t~t + and ti (if kinematically allowed), while direct annihila- 
tion into electron-positron pairs will be strongly suppressed. 
Then we need to consider the following annihilation modes: 
XX -» W + W~, XX -> ZZ, XX -> bb, XX -» t + t~,xX -> it- 
Both the gauge bosons and the fermion pairs produced in 
the neutralino annihilation will initiate a cascade that will 
eventually lead to a continuum of photons, neutrinos, elec- 
tron/positrons pairs and protons final states, extending to 
energies much smaller than the rest mass of the DM particle. 
In Figs. [8] and [9] we show the photon and electron spectra 
dN/ dE for different annihilation channels and for neutralino 
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Figure 8. The photon spectrum from ncutralino annihilation for 
different channels. The solid black lines represent the single anni- 
hilation spectrum for the channel XX ~ * f° r neutralino masses 
of 50, 150 and 600 GeV from left to right. The blue short-dashed 
lines are the same for the \\ — + t+t~ channel. The red long 
dashed lines are for the \\ — > W + W~ channel, for a neutralino 
mass of 150 and 600 GeV. 
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Figure 9. The electron spectrum from neutralino annihilation 
for different channels. The lines are the same as in Fig. \E\ 

masses of 50, 150 and 600 GeV. The sp ectra were calculated 
using PYTHIA0 (|Siostrand et al.ll200ll l. 



The actual spectrum produced by the annihilations will 
depend on the branching ratios of the various channels; this 
in turn will be determined by the gaugino and higgsino frac- 
tions of the neutralino. In the following, we will consider four 
representative supe rsymmetry scenarios , in a similar way to 
what was done by Irlooper et al.l (|2004l ) . First we consider 
a 50 GeV neutralino with an annihilation branching ratio 
of 0.96 to bb and of 0.04 to t + t~ (designated as model 1). 
Such a particle could be gaugino-like or higgsino-like, since 

9 |http:/ /home. thep.lu.se~ torbjorn/Pythia. html 



Figure 10. Photon energy spectrum EdNj/dE from neutralino 
annihilation, for the four models described in the text. 
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Figure 11. Electron energy spectrum EdN e /dE from neutralino 
annihilation, for the four models described in the text. 

for masses below the gauge boson masses, these modes dom- 
inate for either case. 

Second, we consider two cases for a 150 GeV neutralino. 
One (designated as model 2) which annihilates as described 
in model 1 , and another (designated as model 3) which an- 
nihilates entirely to gauge bosons, W + W~ or ZZ . Such 
neutralinos are typically gaugino-like and higgsino-like re- 
spectively. 

Finally we consider heavy, 600 GeV neutralinos, which 
annihilate to bb with a ratio of 0.87 and to t + t~ or t + t~~ 
the remaining time (designated as model 4). 

In Figs. [10] and [TJJ we show spectrum of photons and 
electrons produced in a single annihilation, for these four 
neutralino models. 

Although these models do not fully encompass the ex- 
tensive parameter space available to neutralinos at present, 
they do describe effective MSSM benchmarks. Furthermore, 
the relevant results for neutralinos with a mixture of the 
properties of those above can be inferred by interpolating 
between those presented. 
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We also consider a light dark matter candidate, annihi- 
lating directly into electron-positron pairs. This will result in 
a monochromatic spectrum with E = mnyi- We consider two 
different LDM candidates with masses moM = 3, 20MeV 
respectively. 

The various final annihilation products have different 
interactions with the particles in the IGM. In the following 
we will examine the interaction of photons and pairs with the 
IGM. Protons are very penetrating and thus do not transfer 
energy to the IGM; neutrinos interactions are so weak that 
they are also unable to transfer energy to the IGM, so that 
the annihilation energy that ends up into protons and neu- 
trinos is effectively lost for the purpose of heating/ionising 
the IGM. 



5.1 Photons 

In this paper we are mostly interested in the absorption 
of 7-ray photons. The absorption processes of x-ray and 
7-ray photons at cosmologic al distances were discussed by 
IZdziarski fc Svenssonl (|l989l ). In principle, the possible en- 
ergy loss mechanisms for photons are: photoionisation of 
atoms; Compton scattering on electrons; pair production on 
atoms; pair production of free electrons or nuclei; scattering 
on CMB photons; pair production on CMB photons. The 
total rate for fractional energy loss (j>-f(z, E), i.e., the frac- 
tion of the photon energy that is lost in a unit time, is given 
by a sum over the contributions of the single processes: 



1 dE 
~E~dl 



(41) 



where the index i runs over the different processes listed 
above. However, for z < 500, the dominant processes are 
photoionisation, effective for energies below ~ 10 keV, and 
pair production on CMB photons, effective for energies 
above ~ 10 -f- 100 GeV. Compton scattering is effective only 
for z > 100, in a range of photon energies roughly centered 
around ~ 1 MeV; at z = 500, the region where absorption is 
dominated by Compton scattering extends roughly from 10 
keV to 30 MeV. The other processes, namely scattering on 
CMB photons and pair production on atoms, can be safely 
neglected for our purposes since they are only relevant at 
large redshifts. 

The rate for fractional energy loss by photoionisation 



is given by: 



fa, im (z,E) = T_ v J n b (z)c, (42) 

16 

where rib(z) is the number density of baryons at redshift z, 
and o"He+H is the absorption cross section per helium atom 
(hence the factor of 16 in Eg. 1421 since nn a = n(,/16), given 
by: 



<7He+ H (£0 = 5.1 X 10" 



( E 



V250eV 



(43) 



where p = 3.3 for E > 250 eV, p = 2.65 for 25 < E < 250 eV. 
The fractional energy loss rate by Compton scattering 

is: 



where or is the Thomson cross section, e = E/m e c 2 is the 
photon energy in units of the electron mass, n e ~ 0.88rif, is 
the total electron density at redshift z (including both free 
and bound electrons), and g(e) is: 



9(e) 



(e-3)(e+l) 



ln(l + 2e) 



2 (3 + 17e + 31e 2 + 17e 3 - 10e 4 /3) 
e 3 (l + 2e) 3 



(45) 



The pair production term is: 

7 , P air(2, E) = 27r 1/ V^e 3 /(d9) (46) 

where a is the fine structure constant, A c is the Compton 
wavelength, = kT cm b/m e c 2 is the CMB temperature in 
units of the electron mass, and the function f(x), in the 
limit x <C 1, is given by: 



-1/2 



exp 



l + ?i 
4 



(47) 



^{z,E) = a T eg(e)n e (z)c, 



(44) 



The condition t9 < 1 can be expressed as E <C 10 eV/(l + 
z), so this form of the cross section is adequate for our pur- 
poses. 

In order to assess the efficiency of the energy loss mech- 
anisms, the rate 7 should be compared with the expan- 
sion rate, as given by the value of the Hubble constant 
H(z).When (j>~t 3> H(z), the photon loses all of its energy 
on a time scale small compared to the cosmological time, 
so that the energy loss mechanisms are very effective and 
the Universe is opaque to its propagation. It can then be 
assumed that all the photon energy is instantly deposited in 
the IGM. In the opposite regime, <jf> 7 <C H(z), the photon 
loses a significant fraction of its energy on a time scale larger 
than the cosmological time, and the Universe is effectively 
transparent to the photon propagation. 

This is illustrated in Fig. [12] where we show the pho- 
ton transparency window in the (E, z) plane. For illustra- 
tive purposes, we consider redshifts as large as z = 1000 
in the figure, so that, in addition to the three processes for 
which we have listed explicitly the energy loss rates, we have 
also included the scattering over CMB photons and the pair 
production over nuclei in our computation of the total rate 
7 . In the filled region, 7 > _f/(z), corresponding to the 
optically thick regime. In the white region, </> 7 < H(z), cor- 
responding to the optically thin regime. 



We will compute the absorbed fraction in detail in 
Sec. 15.31 however we can already gain some qualitative in- 
sight by looking at the figure. For the supersymmetric mod- 
els considered, the average energy of the photons produced 
in the annihilation is of order of few GeVs, and in any case 
their energy is at most a few hundreds of GeVs (in the case of 
our heaviest candidate, the 600 GeV neutralino of model 4, 
only ~ 1% of the total energy produced in the annihilation is 
released in the form of photons with energy _E 7 > 200 GeV). 
As it can be seen from the figure, this photons lie in the mid- 
dle of the transparency window: their energy is too low for 
pair production, as already noticed, but on the other hand 
it is too high for photoionisation (or Compton scattering at 
z > 100) to be effective. These photons will propagate freely 
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Figure 12. Photon transparency window. In the black region, the 
photon loses all of its energy through interaction with particle in 
the IGM and CMB photons. In the white region, the photon can 
propagate freely. The dashed lines represent photon trajectories. 



and their energy will decrease due to the cosmological ex- 
pansion. Although it is in principle possible that due to the 
cosmological redshift effect a photon produced in the trans- 
parency window at a given time will be absorbed later, we 
see from the figure that this is practically never the case. 
In conclusion we expect that the absorbed fraction for pho- 
ton will be very small, and that the photons produced in 
neutralino annihilations will instead show up in the diffuse 
gamma background. 

5.2 Electron-positron pairs 

Electrons and positrons can lose energy through collisional 
ionisation of atoms or through inverse Compton scattering 
on CMB photons. In addition, positrons can annihilate with 
thermal electrons, but we have verified that this process is 
always subdominant with respect to inverse Compton scat- 
tering, and we will neglect it in the following. Other energy 
loss mechanisms, like synchrotron radiation loss, can also be 
safely neglected. 

The rate for energy loss through collisional ionisation 

is: 



<pe,ion(E, Z) = 



v 2ire 



E m e v 2 



In 



2T 2 

m e t; 2 7 2 r m ax,Hc 
2J 2 



+ ©(7) 



(48) 



where v is the electron velocity, 7 = E/m e c 2 is the electron 
Lorentz factor, In — 13.59 eV and Thc = 24.6 eV are the 
hydrogen and helium ionisation thresholds respectively, Zu 
and Zhb are the hydrogen and helium atomic numbers, the 
function X>(7) is given by: 

2 



°(7) = TS 



7 T 



I- 1 - 

7 



(49) 



and T maXi H and T maX) He are the maximum energy transfers 
in a single collision: 



T, 



max,H 



27 m H m e u 



m 2 + mfj + 2^/m e mu 
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(50) 
(51) 



The fractional energy loss rate through inverse Comp- 
ton scattering is given by: 



/ „s _ 4 a T d7cMB(.z) 7 2 - 
3 m e c A 7 



1 



(52) 



where (7cmb(z) is the CMB energy density at redshift z. 

In the case of inverse Comptons losses, it must be taken 
into account that the electrons and positrons do not transfer 
their energy directly into the IGM; instead, they accelerate 
the CMB photons, boosting their energy by a factor ~ 7 2 . 
These up-scattered photons can either be absorbed by the 
IGM or escape, depending on their energy. As explained 
above, at redshifts below a few hundreds, the only relevant 
photon absorption processes are photoionisation, compton 
scattering and pair production. The latter is only efficient 
for photon energies below above (at least) ~ 10 GeV; a 
simple calculation shows that the electrons and positrons 
produced in the annihilation of the dark matter candidates 
considered here are not energetic enough to boost the CMB 
photons above the threshold for pair production. We can also 
safely neglect compton scattering, since it is only efficient for 
z > 100 and in a small energy region around 1 MeV. Then 
we need to consider only photoionisation as the secondary 
process leading to the absorption of the photons produced 
by inverse compton scattering of electrons and positrons. 

In order to modelize this effect, we proceed as fol- 
lows. We consider that a scattered photon is absorbed 
through photoionisation if 7 ,i O n(z, E) > H(z); this con- 
dition amounts to the requirement that the energy of the 
photon after the interaction should be smaller than 



(z) = 0.64 keVx 



0.11 



VL b h 2 
0.0224 



(1 + 2 



,3/2 



(53) 

where q — 0.3, and we have assumed a matter dominated 
Universe. 

The energy density f/rjMB °f photons that will be ab- 
sorbed after the interaction, as a function of the electron 
Lorentz factor and of redshift, is then given by: 



^cmb(7, 2 ) = 



8tt 
(he) 3 



M o»/7 2 



dE- 



E 3 



io„/7 2 



£/fcT CMB M - 1 : 



(54) 

where the factors of 7 2 in the integration limits are due to 
the fact that the energy of the photon is increased by a 
factor of 7 2 after the interaction. The above expression can 
be recast as follows: 



U'c 



CMB 



UCMB f 
7T 4 /15 



e» - 1 



dy, 



(55) 



where y mils = B m in,ion/(7 2 fcTcMB), and analogously for 

Finally, we can estimate the rate of energy injection into 
the IGM by inverse Compton $ eiCom (z, E) by replacing the 
total energy density Ucmb in Eq. (|52p with the "reduced" 
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Figure 13. Transparency window for electrons. In the white re- 
gion, electrons propagate freely. In the gray regions, electrons 
transfer all of their energy to the CMB photons, but these are 
subsequently lost, so that no energy is injected in the IGM. In 
the black regions, all the electron energy is efficiently injected in 
the IGM. See text for discussion. 



energy density U' CMS . This amounts to weighting the inverse 
compton absorption rate <jie,com by an efficiency factor 77 = 
^cmb/^CMb: 



l(z,E) = T](z,E) </> e ,com(z, E). 



(56) 



The behaviour of electron-positron pairs with respect 
to the energy transfer to the IGM is summarized in Fig. 
1131 In the white region, the total energy loss rate is smaller 
than the expansion rate: e ,ion + </>e,com < H , so that Uni- 
verse is transparent to the propagation of electrons. In the 
grey regions, the electrons and positrons interact by inverse 
Compton scattering, but the resulting photons fall in the 
photon transparency window. This means that the Uni- 
verse is opaque to the propagation of electrons, but nev- 
ertheless their energy is not transferred to the IGM. This 
corresponds to the regime in which <f> e ,com > H, ^ e ,ion but 
$ e , com < H. Finally, in the black regions the electron en- 
ergy is efficiently transferred to the intergalactic medium. 
In particular, the region on the left correspond to the case 
in which e ,ion > H, e ,com, meaning that collisional ionisa- 
tion is the dominant process; the region on the right corre- 



spond to 



and <E> e 



> H, meaning that 



inverse Compton is the dominant interaction, and that the 
upscattered CMB photons fall into the photon absorption 
window. 



5.3 Energy absorbed fraction 

The energy deposition of p hotons and pairs in the IGM 
has be en worked in detail bv lRipamonti. Mapelli fe Ferraral 
(|2007l ). Here we will consider a more simplified approach, 
in the form of a "on the spot" approximation. In practice, 
we will be comparing the energy loss rate <j> of the product 
particle, for example a photon, with the cosmological expan- 
sion rate H. If <f> > H, we consider that all the energy of the 
particle is immediately absorbed by the IGM; if on the con- 
trary <f> < H, we consider that the particle energy is "lost" 



for the purpose of heating the IGM. In the case of electrons 
and positrons losing energy via inverse Compton scattering 
with the CMB photons, we take into account the efficiency 
related to the transfer of energy from the scattered photons 
to the IGM, as explained in the previous section. 

In more detail, what we do is the following. At every 
redshift, starting from z = 500, we compute the absorp- 
tion region Ai(z) for every kind of particle, i.e., the energy 
range in which particles of that kind efficiently transfer their 
energy to the IGM, using the (f> «s H criterium. We note 
that the region A% will possibly consist of two or more dis- 
connected intervals. For example, by looking at Fig. 1131 we 
see that for electrons the region A e (z = 100) is roughly 
the union of the two intervals 10 eV $C Ek ^ 100 keV and 
lOMeV < Ek 1 GeV, Ek being the kinetic energy of the 
electron. 

Then we evaluate the energy absorbed fraction simply 
as the ratio of the energy that is absorbed to the total energy 
produced in the annihilation,: 



E 



/abs (2) 



E d Jm dE 

dE 



(57) 



2m dm c 2 

where the sum runs over the different annihilation prod- 
ucts. For electrons and positrons, in the region where inverse 
compton scattering is the dominant loss mechanism, we ad- 
ditionaly weight the integrand in Eq. (I57|l with the efficiency 
factor 77 defined in the previous section, in order to account 
for the finite fraction of the energy lost by the electron that 
is actually injected into the IGM. 

In Fig. [14] we show the energy absorbed fraction for the 
supersymmetric models described at the beginning of this 
section. We see that for all four models, / a b s is nearly con- 
stant and always much smaller than unity at the redshifts 
of interests. In particular, / a b s goes from ~ 10 -2 for the 
model with the lightest candidate (model 1, m x = 50 GeV) 
to ~ 10~ 4 for the model with the heaviest candidate (model 
4, m x — 600 GeV). This behaviour can be understood as 
follows. The photons and electrons produced in neutralino 
annihilations peak at E ~ 1 GeV for models 1 to 3, and at 
E ~ 10, GeV for model 4. Thus photons are in the middle 
of the transparency window and do not get absorbed. Elec- 
trons (and positrons), on the other hand, possibly lie in (or 
close to) the region in which energy is transferred to the IGM 
by means of inverse Compton scattering over CMB photons 
and the subsequent absorption of those through photoionisa- 
tion. This region is roughly centered around E ~ 0.1 GeV so 
that lighter neutralinos, producing softer annihilation spec- 
tra, will be absorbed more easily. This also explains while 
model 3 neutralinos inject more energy than their model 2 
counterparts: the neutralino mass is the same in both mod- 
els, both the former produces sligthly less energetic elec- 
trons. 



In Fig. [15] we show the energy absorbed fraction for the 
two LDM models considered. In this case we see that the 
energy absorbed fraction is larger with respect to the SUSY 
scenarios. For the lightest candidate, ttildm = 3, GeV, the 
energy absorbed fraction is 1 at z = 1000, then decreases 
quite rapidly, reaching a minimum value / a b s — 5 x 10~ 2 at 
redshift 100, and again starts to increase slowly. In the case 
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Figure 14. Energy absorbed fraction / a b s for the four neutralino 
models described in the text. The red solid line is for a 50 GeV 
neutralino. The green long-dashed line is for a 150 GeV, gaugino- 
like neutralino. The blue short dashed line is for a 150 GeV, 
higgsino-likc neutralino. The magenta dotted line is for a 600 
GeV, gaugino-likc neutralino. A smaller neutralino mass leads to 
a larger absorption. See text for discussion. 

of the heaviest candidate, ttildm = 20 GeV, the absorbed 
fraction stays constant and equal to 1 until z ~ 10, where 
it starts to decrease. This behaviour can also be understood 
in terms of the transparency/absorption windows studied in 
the previous section. We considered LDM decaying directly 
to electron-positron pairs only, so that every particle pro- 
duced carries an energy equal to the DM particle mass. As it 
can be seen from Fig. 1131 a 3MeV electron lies in the absorp- 
tion region only if it is produced at redshifts close to 1000. 
Electrons produced later lie outside this region, and this ex- 
plains why / a bs decreases between z — 1000 and z = 100. 
The increase after z = 100 is due to the increased impor- 
tance of collisional ionisation (the ratio (f>e,ioii/H scales with 
redshift like (1 + zf /2 ). On the other hand, 20 MeV elec- 
trons always lie in the absorption region from z = 1000 to 
z — 10, and only after that they fall ouside it; this is why 
the absorbed fraction is always ~ 1 in this redshift range. 



In the derivation above we have implicitly assumed 
that the Universe is mostly neutral in the redshift range 
of interest. When the Universe starts to be reionised, the 
energy loss rate for photoionisation will become smaller, 
since it is proportional to the density of neutral atoms. 
From the discussion above, we have seen that the energy 
produced in DM annihilations is mainly injected into the 
IGM through inverse Compton scattering over CMB pho- 
tons (e7cMB — > &i) followed by absorption of the photon 
via photoionisation (jA — > A + e~). Then we expect the ab- 
sorbed fraction to start going rapidly to zero once reionisa- 
tion begins, at z ~ 10. 



6 THE 21 CM BACKGROUND 

The emission or absorption of the 21cm line signal ema- 
nating from neutral gas is associated with the transition 



Figure 15. Energy absorbed fraction / a b s for the two LDM mod- 
els described in the text. The red solid line is for a 3 MeV LDM 
particle. The green long-dashed line is for a 20 MeV LDM parti- 
cle. A larger particle mass leads to a larger absorption. See text 
for discussion. 

between the n — 1 triplet and singlet hyperfine levels of 
hydrogen. The transition rate is governed by the spin tem- 
perature, T a , defined as 

^=3exp(^V (58) 
no \ Is J 

where no and rii are the respective number densities of hy- 
drogen atoms in the singlet and triplet states, and T* = 
0.068K is the equivalent temperature corresponding to the 
transition energy. 

6.1 CMB - kinetic temperature coupling 

In the presence of the CMB radiation field, the spin tem- 
perature of the neutral hydrogen gas rapidly tends towards 
the CMB temperature T C mb ^ 2.725(1 + z)K. In order for 
neutral hydrogen gas to produce a detectable signal in the 
21cm background, be it absorption or emission, that is dis- 
tinguishable from that generated from the CMB, the kinetic 
temperature Tk of the gas must decouple from Tomb • 

In a Universe containing stable, non-annihilating dark 
matter, the spin temperature and the kinetic temperature 
of HI gas are coupled to the Tcmb until z ~ 200. At 
30 < z < 200, prior to the formation of non-linear baryonic 
structures, the IGM cools adiabatically, i.e. Tk oc (1 + z) 2 , 
compared to Tcmb oc (1 + z). During this epoch spin- 
exchange collisions between hydrogen atoms, protons and 
electrons are efficient at coupling Tk and Ts of the gas, and 
consequently an absorption at wavelength A = 21(1 + z) cm 
can be observed until approximately z ~ 70. At later times 
cosmological expansion reduces the frequency of these col- 
lisions significantly to the extent where Ts re-couples with 
Tcmb, diminishing the 21cm absorption signal. 

However, in a Universe containing annihilat- 
ing/decaying dark matter which injects an appreciable 
of energy into the IGM, the thermal history of the gas 
may be significantly altered to the extent where the 
corresponding changes in the evolution of the 21 cm signal 
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are detectable by current and future radio interferometers. 
Of particular importance is the high sensitivity of these 
changes with respect to the nature of the dark matter, 
making them a powerful tool for constraining the properties 
of potential dark matter candidates. 

There are two mechanisms which can decouple Ts from 
Tomb: firstly, the fore-mentioned spin-exchange collisions 
between neutral atoms, electrons and protons, which are ef- 
fective at z ^ 70 before the Hubble expansion has rarefied 
the gas in the IGM, and secondly, scattering by Lyman-a ra- 
diation ( "Wo uthuysen-Field" effect, also known as "Lyma n- 
q pumping" ijWouthuvsenll 19521 ; iFieldl Il959l : lHiratall2006h L 
which couples Ts to Tk via the mixing of the n = 1 hyper- 
fine states through intermediate transitions to the 2p state. 

In the quasi-static approximation for the population of 
the hyperfine levels in question, and in the absence of ra- 
dio sources, the HI spin temperature is a weighted mean 
involving T K and Tomb, 



Tcmb + vTk 



1 + y 

The coupling coefficient y can be written as 

y = y a + yc, 



(59) 



(60) 



where y a is the term associated with Lyman-a pumping, 
given by 

y a = , (61) 

AiqI K 

whilst yc is associated with the de-excitation of the HI hy- 
perfine triplet state due to collisions with neutral atoms, 
electrons and protons, collectively written as 



6.2 Modifications to IGM thermal evolution when 
incorporating DM 

In this section we describe the modifications to the stan- 
dard equations describing thermal and ionisation history of 
the IGM when we incorporate the potentially significant en- 
ergy deposition of the products of annihilating dark matter, 
which in turn determines the 21 cm signal. 

We parameterize the effect of dark matter annihilation 
by the rate of energy injection given by eq. (JTJ. This en- 
ergy is then used to excite and ionise the hydrogen and 
helium atoms in the IGM. We will not enter in the de- 
tail of the partition of energy between hydrogen and he- 
lium, but instead assume that it is divided proportionally 
to the respective number densities. This means that a frac- 
tion 1/(1 + /He) of the absorbed energy will go to hydrogen, 
while a fraction /h c /(1 + /h c ) will go to helium, /He be- 
ing the helium to hydrogen number ratio. Then we need to 
know how the energy is partitioned between the different 
processes. The relative fractions x»>Xh an d Xe of the en- 
ergy absorbed which is diverted towards ionising, heating 
and exciting hydrogen and he lium atoms were calculated by 
IShull fc van Steenberd (Il985l). Their resu lts can be approx- 
imated by ( Chen fe Kamionkowskil|2004l 'l 



XiA z ) 
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(64) 
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where Xj (z) is the ionisation fraction of the relevant nuclear 
species j (i.e. j=H or He for hydrogen or helium nuclei re- 
spectively), defined as 



yc 



A w Tk 



(C H + C e + C p ). 



(62) 



In the above equations Aio = 2.85 x 10~ 15 s is the rate 
of spontaneous photon emission, Piq is the de-excitation 
rate of the hyperfine triplet state due to Lyman-a scat- 
tering and Cn,C e and C p are the de-excitation rates as- 
sociated with collisions of hydrogen atoms with other hy- 
drogen atoms, electrons and protons respectively. We write 
Pio = {lQ-Kj a cJa) I '(27h p v a ), where J a is the background 
intensity of Lyman-a photons, a a is the Lyman-a photon 
absorption cross-section for neutral hydrogen and h p = 6 is 
Planck's constant. We neglect the small co rrections to the 
above expressions proposed bv lHiratal (|2006T )). The H-H col- 
lision term can be written as Ch = fcionm, where kio is the 
effective single-atom collision rate coefficient for which we 
adopt the fit fcip = 3.1 x 10~ 11 r^ 357 exp(-32/r A -) cimV 1 
proposed by iKuhlen et al.l (|2005l ). which is accurate to 
within 0.5% in the range 10 < T K < 10 3 K. For the e- 
H collision ter m, C e = n e 7 e, we have used the following 
fit proposed by iLisztl (|200ll ). log^e/cmV 1 ) = -9.607 + 
0.51og(T x )exp [-(logT K ) 4 - 5 /1800] for 1 < T K < 10 4 K 
log(7 e /cm 3 s~ 1 ) = -9.607 + 0.51og(T K -) for T K < IK (|Smithl 
1 19661 ). and j e (T K > 10 4 K) = 7e (10 4 K). We ignore de- 
excitations involving collisions with protons since they are 
typically much weaker than those involving electrons at the 
same temperature. 



(66) 



where is the number density of ionised nuclei of the 
species j. We can also define a total ionisation efficiency 
Xt = (Xi,n + /hcX»,Hc)/(1 + /He), and similar quantities for 
heating and excitation. 

We compute ionisation and thermal history of the 
IGM when incorporating our chosen species of dark 
matter using the publicly availab le code RECFAST 
|Seager. Sasselov fc Scott) 1 19991 . l200Ch . modifying the stan- 
dard evolution equations for the ionisation fractions 
of hydrogen and helium nuclei, aswell as the evolu- 
tion equation for the kinetic t emperature as follows 
(|Padmanabhan fc Finkbeinerll2005l '). 
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A further equation needed to calculate the 21 cm sig- 
nal is that describing the evolution of the Lyman-a back- 
ground intensity J a , which can couple the spin and ki- 
netic temperatures of the H-atoms in the IGM through the 
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Wouthuysen-Field effect. H-atoms, excited by collisions with 
fast photoelectrons subsequently produce a cascade of line 
photons, including Lyman-a photons which are then likely 
to be re-absorbed by the opt ically-thick IGM . We u tilise 
the approximation adopted by iFurlanetto et alj (|2006f ) that 
approximately half of the total energy which is diverted to 
excite hydrogen is used to produce Lyman-a photons, i.e. 

Xa ~ Xe,H/2 ' I ; 1 J n 

Following the treatment bv lValdes et al.l (|2007l ) we ob- 
tain (however, please note the additional factor of v a in front 
of the expression) 



n\hcv a 
4nH(z) 
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X e X e ,na 2 2 p + £e,H:rHl7cH H ; 

nuhv a 



(70) 



where the first two terms are the contributions are as- 
sociated with the collisional excitation involving electrons 
discussed above, and the last term is the contribution 
from dark matter. Also, ct^fp is t he effective rec ombi- 
nation coefficient to the 2 2 P level |Pengellvl [l964h , and 
7 eH ~ 2.2 x 10 _8 exp [-11.84/(r/10 4 K)]cm 3 s _1 is the col- 
lisional excitation rate of H I atoms involving electrons 
jShull fc van Steenberdll985h . 

The quantity most intimately associated with the obser- 
vations of the cosmological 21 cm signal intensity is the dif- 
ferential brightness temperature deviation, 8Tb, betwe en the 
21 cm signal and the CMB, appro ximately given by (|Fieldl 
1 19581 . 1 19591 ; ICiardi fc Madaull2003l ) 
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(71) 



where x^i = 1 — £h is the averaged fraction of neutral hy- 
drogen in the patch of sky being observed. 



7 RESULTS 

In the following section we illustrate our predictions for the 
effects on the thermal history of the IGM caused by the 
additional energy injected into it by annihilating neutralino 
and light dark matter, when including the effect of the en- 
hancement from dark matter structures. 

Since, as we have seen, this enhancement can be very 
large, boosting the DM annihilation rate by several orders 
of magnitude, we want to be sure that this does not spoil 
other observations. Then we perform two tests on each of the 
clumping factors investigated, before taking it into consider- 
ation for the purpose of its effect on 21cm observations. First 
of all, we check that the huge injection of energy into the 
IGM does not lead to a early ionisation. We use to this pur- 
pose our modified version of the RECFAST code described 
before, and discard all the clumping factors for which the 
ionised fraction x e > 0.01 at z = 14. We choose this value 
of the redshift because it is close to the 2a upper limit 
to ^re ion coming from WMAP5 observations (|Dunklev et al.l 
120081 ). Secondly, we check that the diffuse photon flux pro- 
duced does not exceed observed diffuse background (adopt- 
ing the conservative approximation that / a b s ~ 0, that for 
all the models we consider is quite a good approximation at 



the present time z = 0, where the clumping factor reaches 
its maximum value). We use to this purpose the measure- 
ments of the diffuse gamma-ray backgrou nd in the 1 MeV 
- 100 GeV range conducted by EGRET (|Sreekumar et al.1 
ll99StlStrong. Moskalenko fc Reimerll2004l ) and COMPTEL 
|Weidenspointner et al.ll2000l ). and those of the diffuse X- 
ray background in the sub-MeV ra nge conducted by the SPI 
spectrometer aboard INTEGRAL (|Churazov et aLlBoOTT ). 

In the following, owing to the fact that the effects on 
the spin and brightness temperature can be very subtle, we 
display results only for the the most optimistic clumping 
factors that conform to the above criteria. For reference, in 
the appendix to this paper we have tabulated the relevant 
astrophysical parameters associated with all clumping fac- 
tors investigated, indicating which have been excluded on 
the basis of the criteria described above. 



7.1 Neutralino dark matter 

We begin with an analysis of neutralino dark matter. Whilst 
in all four of the models considered the deviations of the 
spin temperature from that calculated in a scenario where 
dark matter is excluded is marginal, there are significant 
deviations in the evolution of the kinetic temperature which 
we highlight below. 

We start with model 1, a 50 GeV neutralino with a 
substantial Higgsino and Gaugino fraction that annihilates 
to bb pairs 96% of the time and to t + t~ otherwise, with 
a canonical annihilation cross-section of (<7 ann .i;) = 3 x 
10~ 27 /f2 D M,ofr 2 2± 2.7 x 10~ 26 cm 3 s _1 . Owing in particular to 
the very light nature of this neutralino, the energy injection 
rate is very large (since overall, edm scales as w,q M ). Con- 
sequently, the majority of the clumping factors calculated 
using the Moore profile are excluded owing to the overpro- 
duction of gamma-rays. In figure[TS] we display the effects 
of such dark matter particles on the evolution of the Tk 
and Ts (dot-dashed and dashed curves respectively) for the 
most optimistic clumping factors consistent with observa- 
tions of the diffuse background, namely M18 (blue) and N4 
(magenta) calculated using Moore and NFW density pro- 
files respectively, compared to a scenario where dark matter 
is absent (black dot-dashed and black solid curves respec- 
tively). As can be seen from figure [T6l significant heating of 
the IGM by dark matter begins when the respective clump- 
ing factors start to elevate significantly above unity, which 
for M18 corresponds to z — 20 and z — 90 for N4, indicat- 
ing approximately when the lightest dark matter structures 
start to form. 

Next we consider dark matter composed of neutralinos 
described by model 2, that is, 150 GeV Gaugino-dominated 
neutralinos, that annihilate 96% to bb and 4% to t + t~. Once 
again, the majority of the clumping factors calculated using 
the Moore profile are excluded owing to the overproduc- 
tion of gamma-rays, with M8 being the most optimistic non- 
excluded clumping factor, whilst for structures with NFW 
profiles all but Nl, N5 and N9 are permitted, with N2 be- 
ing the most optimistic. In figure fTTl we display the effects 
of such dark matter particles on the evolution of the Tk 
and Ts (dot-dashed and dashed curves respectively) for the 
clumping factors M8 (blue) and N2 (magenta) calculated us- 
ing Moore and NFW density profiles respectively. From fig- 
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Figure 16. The evolution of the IGM kinetic and spin tempera- 
ture (dot-dashed and dashed curves respectively) for the clump- 
ing factors M18 (blue) and N4 (magenta) calculated using Moore 
and NFW density profiles respectively, for structures composed 
of 50 GeV ncutralino dark matter described by model 1 with an- 
nihilation cross-section (<T ann .u) = 2.7 X 10 — 26 cm 3 s _1 , compared 
to a scenario with no dark matter (black dot-dashed and black 
dashed curves respectively). Also displayed is the temperature of 
the CMB (red solid curve). 
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Figure 17. Same as for figurc lTSI but for dark matter composed 
of 150 GeV Gaugino-dominated ncutralinos described by model 2 
and using clumping factors M8 (blue) and N2 (magenta). 



ure !17l the additional heating of the IGM by dark matter be- 
gins at similar times for both profiles, at around z — 40, al- 
though both M18 and N2 start to elevate significantly above 
unity at around z = 100. At all except recent times M8 is 
greater than N2 explaining the consistently increasing rate 
of T K for M8 relative to that for N2. 

Next we consider dark matter composed of neutralinos 
described by model 3, that is, 150 GeV Higgsino-dominated 
neutralinos, that annihilate 58% to W + W~ and 42% to 
ZZ. One would expect that since the results associated with 
such dark matter to be very similar to that for neutralino 
dark matter described by model 2 since the energy injec- 
tion rate is largely set by the dark matter particle mass 
m x . However, clumping factors M7 and M14 for structures 
possessing Moore profiles and N5 for NFW profiles are now 
permitted, because the annihilation products resulting from 
model 3 neutralinos produce an overall lower flux in pho- 
tons. In figure lrSl we display the effects of model 3 neutralino 
dark matter on the evolution of the Tk and T$ (dot-dashed 
and dashed curves respectively) for the clumping factors M7 
(blue) and N2 (magenta). The similarity of these effects com- 
pared to the 150 GeV Gaugino-dominated neutralinos de- 
scribed by model 2 reaffirms the notion that it is the largely 
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Figure 18. Same as for figurc lTEl but for dark matter composed 
of 150 GeV Higgsino-dominated neutralinos described by model 3 
and using clumping factors M7 (blue) and N2 (magenta). 




Figure 19. Same as for figurc lTEl but for dark matter composed 
of 600 GeV Gaugino-dominated ncutralinos described by model 4 
and using clumping factors M4 (blue) and N2 (magenta). 



the neutralino mass that dictates the resulting deviations 
in the 21 cm signal rather than its Gaugino fraction. This 
result is not altogether surprising if one recognizes that the 
majority of the energy injected into the IGM by neutralinos 
ultimately ends up as electrons, positrons and photons re- 
gardless of the initial annihilation modes to the extent where 
these various subtleties between different types of neutralino 
are largely lost in the evolution of the Tk and T3. 

Finally, for dark matter composed of neutralinos de- 
scribed by model 4, that is, 600 GeV Gaugino-dominated 
neutralinos that annihilate 87% to bb and 13% to t + t~. 
In figure[Tj|]we display the effects of model 4 neutralino dark 
matter on the evolution of the Tk and Ts (dot-dashed and 
dashed curves respectively) for the clumping factors M4 
(blue) and N2 (magenta). All deviations from the "no dm" 
scenario are extremely suppressed even for the most opti- 
mistic clumping factors, which in this case are M4 and N2, 
for structures possessing Moore and NFW profiles respec- 
tively. This is due to the extremely suppressed energy in- 
jection rates arising from the large mass of the neutralino, 
relative to other neutralinos possessing similar clumping fac- 
tors, with Tk only beginning to deviate from the "no dark 
matter" model at recent times when the clumping factors 
are largest. 

7.2 Light dark matter 

We now turn our attention to LDM. We consider two 
models, firstly dark matter composed of 3 MeV LDM par- 
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Figure 20. The evolution of the IGM kinetic and spin tempera- 
ture (dot-dashed and dashed curves respectively) for 3 McV (blue) 
LDM particles using the clumping factor B12 and an annihilation 
cross-section (crv) = 1.29 X 10~ 29 cm 3 s _1 , and also for 20 MeV 
(magenta) LDM particles using the clumping factor B7 and an an- 
nihilation cross-section (crv) = 8.60 X 10~ 29 cm 3 s~ 1 . The results 
for the evolution of the kinetic and spin temperature are com- 
pared to a those calculated when dark matter is absent (black 
dot-dashed and black dashed curves respectively). Also displayed 
is the temperature of the CMB (solid red curve). 



tides with an annihilation cross-section (crv) = 1.29 x 
10 _29 cm 3 s _1 , and secondly dark matter composed of 
20 MeV LDM particles with an annihilation cross-section of 
(crv) = 8.60 x Hr 29 cm 3 s _1 . 

In figure [20l we display the effects of these dark matter 
candidates on the evolution of the kinetic and spin temper- 
ature of the IGM (dot-dashed and dashed curves respec- 
tively for 3 MeV (blue) and 20 MeV (magenta) LDM parti- 
cles). Once again, for each model we use the most optimistic 
clumping factors consistent with observations of the diffuse 
background. For 3 MeV particles this corresponds to B12, 
calculated using a Burkert density profile. For 20 MeV par- 
ticles the absorbed fraction, fats, of energy injected into the 
IGM is significantly larger than for the lighter candidate, 
allowing for larger enhancements to its rate of annihilation 
by dark matter over-densities, with B7 being the most opti- 
mistic clumping factor permitted. 

As can be seen from figure [20l there is a marked differ- 
ence in the effects of the two models. Both the 3 MeV and 
20 MeV LDM models significantly heat the IGM at early 
times, resulting in Tk decoupling from its standard adia- 
batic revolution following the formation of the first struc- 
tures, which corresponds to 10 6 Mq halos at z ~ 50, and 
46 Mq at z ~ 40 respectively. For 3 MeV LDM this produces 
only modest changes in Ts relative to the standard case, 
slightly reducing the absorption signature approximately in 
the range 50 < z < 100. However for 20 MeV LDM the 
larger clumping factor and absorbed fraction results in in- 
creasingly larger Tk to the extent where Ts exceeds Tomb 
at z ~ 50, resulting in an subsequent emission signature for 
the 21 cm signal that is still significant at z = 10. 



7.3 The 21cm global signature 

Using the above results for the evolution of Ts, calculated 
for our neutralino and LDM candidates, we now calculate 
the evolution of the differential brightness temperature 8Tb 
associated with the cosmological 21 cm signal , which is to be 



directly meas u red b y radio interferometers such as LOFAR 
(see e.g. lBestl (120081 )1. 

In figure l2"T1 we display our results for the predicted evo- 
lution of 8Tb for dark matter composed of annihilating LDM 
particles of mass 3 MeV (blue dashed) and 20 MeV (ma- 
genta dashed) relative to the standard scenario where dark 
matter is absent (black solid). For clarity, we neglected to 
superimpose results for 8Tb associated with neutralino dark 
matter on the same axes owing to the very small devia- 
tions of these results from the standard case 8Tb, o, relative 
to those obtained for LDM. However, in figure [22l we have 
displayed results for the deviations 8Tb — 8Tb,o for neutralino 
dark matter, and again for LDM displayed separately in fig- 
ured 

Firstly considering neutralino dark matter, the evo- 
lution 8Tb is rather similar for all four models, differ- 
ing approximately only in an overall normalisation result- 
ing from competing effects involving a relative enhance- 
ment/suppression owing smaller/larger particle masses and 
larger/smaller clumping factor. As mentioned above, values 
of the brightness temperature deviation \8Tb — 8Tb : o\ relative 
to the standard scenario remain very modest, with maxi- 
mum values of ^ (2 — 3) mK, at the two characteristic peaks 
occurring at z ~ 50—70 and z ~ 20 — 30. However, at z ~ 10, 
such deviations are significantly reduced in all cases to less 
than 0.5 mK. 

For LDM, the situation is far more interesting than for 
neutralino dark matter, since values of the brightness tem- 
perature deviation are far more significant. The evolution of 
8Tb is almost identical for both LDM for z > 100, with devia- 
tion from the standard case increasing up to ~ 15 mK as the 
energy injected into the IGM by LDM annihilations forces 
8Tb to -30 mK at z = 100, compared to -50 mK when dark 
matter is absent. Subsequently, the rapidly decreasing value 
of /abs. for 3 MeV particles results in a steadily decreasing 
value of 1 8Tb — 8Tb,o\, reaching zero at z ~ 30 and despite 
a slight increase in its subsequent evolution, \8Tb — <STi,,o| 
returns to almost zero at z = 10. However, for 20 MeV 
LDM particles, / a b s . maintains its large value until much 
later times resulting in the continued increase in 8Tb — 8Tb,o 
beyond z = 100, reaching a maximum of ~ 35 mK at z ~ 40. 
The subsequent rapid decrease in / a b s . once again results in 
the steady decrease of 8Tb ~ 8Tb ,o, but the rate of decrease 
is sufficiently slow such that 8Tb — 8Tb,o — 9 mK at z — 10, 
which should be comfortably detectable by interferometers 
such as LOFAR despite the various forms of foreground con- 
tamination (see §[8]). However, as already discussed at the 
end of Sec. [5] this result holds if the Universe is still mostly 
neutral at this redshift. If this is not the case, we expect the 
signal associated to DM annihilations to rapidly go to zero 
as the Universe gets reionised. 



8 DISCUSSION 

We have calculated predictions for the effects on the evo- 
lution of the cosmological HI 21 cm signal during the Dark 
Ages for various forms of annihilating neutralino and light 
dark matter. Our calculations differ significantly compared 
to those conducted in previous studies since we include the 
enhancements in the annihilation rate arising from dark 
matter structures. We utilised results from the most recent 
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Figure 21. The 21cm differential brightness temperature as a 
function of redshift. The black solid line displays the evolution 
of 8Tb for the standard scenario where dark matter is absent, 
while the dashed lines correspond to 8Tb for a Universe containing 
LDM particles of masses 3 MeV (blue dashed) and 20 MeV (ma- 
genta dot-dashed). Calculations were performed using the most 
optimistic clumping factors associated with structures possessing 
Burkert profiles that are consistent with the diffuse X-ray back- 
ground, which in this instance were B7 and B12 for 3 MeV and 
20 MeV LDM particles respectively. 




Figure 22. Difference in the 21cm differential brightness tem- 
perature <5Tj, between that calculated for the standard case and 
that calculated using neutralino dark matter described by model 1 
(blue), model 2 (red), model 3 (green) and model 4 (magenta). 
The results displayed correspond to those calculated using the 
most optimistic clumping factors for structures possessing Moore 
(solid) and NFW (dashed) profiles that were used to generate the 
results presented in §[7] 
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state-the-art N-body simulations of dark matter halos to cal- 
culate the evolution of the enhancement in the annihilation 
rate, referred to as the "clumping factor" , owing to the dis- 
tribution of halos and the self-similar distribution of increas- 
ingly smaller substructures predicted to exist within them, 
for a range of values of astrophysical parameters consis- 
tent with the uncertainties in the dynamics of the simulated 
halos. We performed detailed calculations of the absorbed 
fraction of the energy injected into the IGM by the anni- 
hilation products of our dark matter candidates. We used 
the standard equations for the evolution of the kinetic and 
spin temperatures of the IGM with modifications to account 
for the additional energy injected into, and subsequently ab- 
sorbed, by the IGM from dark matter annihilations. Finally, 
we calculated the resulting deviations in the evolution of the 
differential brightness temperature 8Tb, which is the quan- 
tity directly observed by radio interferometer experiments, 
relative to a scenario where dark matter is absent. If such 
deviations are sufficiently large, measurements by future ra- 
dio telescopes, such as MWA, 21 CMA and SKA, may be 
sufficiently sensitive to discern and distinguish the effects of 
different species of dark matter. 

Of particular interest are the impending LOFAR epoch 
of reionization experiments (LOFAR-EoR hereafter) de- 
signed to observe radio fluctuations at frequencies 115- 
215 MHz, corresponding to the reds hifted 21 cm signal 
in the range 6 < z < 12 (see e.g. iBestj ( |2008t )). Un- 
fortunately, the cosmological signal is contaminated by 
a plethora of astrophysical and non-astrophysical com- 
ponents - including, Galactic synchrotron emiss ion from 
diffuse and localized sources (IShaver et al.lll999T l. Galac- 
tic free-free emission (|Shaver et al.lll999T ), integrated emis- 
sion from extragalactic sour c es, such as radio galaxi e s and 



sion from extragalactic sourc es, such as radio galaxi e s and 
clusters (IShaver et al.l Il999l; |Pi Matteo et ahl |2002| . l2004j ; 
lOh fc Mackll2003l ; ICoorav fc Furlanettd I2004T ), ionospheric 
scintillation and instrumental response - the fluctuations of 
which can be up to thre e orders of magn itude greater than 
the cosmological signal (| Jelic et al.ll2008h . 

It is beyond the scope of this paper to provide a thor- 
ough treatment of the above forms of contamination, their 
modeling, and the methods employed to ultimately extract 
the cosmological signal. Many authors have studied these 



the cosmological signal. Many authors have studied these 
topics in detail (IShaver et alii 19991: iDi Matteo et al.ll2002l; 
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Zaldarriaga et all 120041: ISantos et al.l 120051: iMorales et all 
20061 : 1 Wang et a l. 2006; Gleser et al.ll2007l : I Jelic et al.ll200St ) 
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and we refer the reader to such papers for further informa- 
tion. Here, we simply assume that such methods will be 
reliable enough so that LOFAR possesses a sensitivity in 
measurements of 8Tb, after signal extraction, of order 1 mK 
for redshifts z < 12. |^j These estimates are based on prelimi- 
nary results derived from simulations, conducted by LOFAR 
scientists, investiga ting the removal of foregrounds from the 
21cm background (| Jelic et al.ll2008T l. 

We now apply the above sensitivity estimates to the 



Figure 23. Difference in the 21 cm differential brightness tem- 
perature 8Tb between that calculated for the standard case and 
that calculated for 3 MeV (blue dashed) and 20 MeV (magenta 
dot-dashed) LDM as displayed in figure [211 



10 We should mention however that the first luminous astrophys- 
ical objects in the Universe forming near to the epoch of reioni- 
sation, which from the most recent measur ements by WMAP in - 
dicate could occur as early as z = 14 (2 a) llDunklev et alj '2008). 
could also produce a significant contribution to the global 21cm 
signal. 
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results from our analysis. Firstly, we investigated dark mat- 
ter composed of neutralinos, possessing masses ranging from 
50-600 GeV, with being some Gaugino/Higgsino dominated, 
and an annihilation cross-section (av) ~ 2.7 x 10 _26 cm 3 s _1 
which is consistent with current estimates for the currently- 
inferred abundance of dark matter in the Universe. Our re- 
sults indicate that the evolution of the brightness temper- 
ature deviations 5Tb — 5Tb, a for the widely varying forms 
of neutralino dark matter considered are similar to within 
a normalisation factor of order unity. Deviations of up to 
±2 — 3 mK are predicted to occur in the range z — 50 — 70 
and z = 20 — 30. However, at frequencies detectable by LO- 
FAR, such deviations are significantly reduced in all cases 
to less than 0.5 mK at z ~ 12, placing them on the very 
edge of LOFARs detectability threshold, but certainly may 
be detectable by future experiments. 

We then investigated the effects of dark matter com- 
posed of annihilating LDM with masses of 3 and 20 MeV, 
with annihilation cross-sections of (av) ~ 1.29 x 10 -29 and 
8.60 x 10 -29 cm 3 s _1 , based on constraints derived from 
their predicted effects on the CMB. Our results predict sig- 
nificant deviation for 3 MeV LDM particles in the range 
15 < z < 600, rising to as much as +15 mK at z ~ 100. How- 
ever at later times, owing to the rapidly decreasing rate at 
which the IGM absorbs energy from the annihilation prod- 
ucts of these particles, deviations in the brightness temper- 
ature are almost negligible at z = 12 and are unlikely to 
be detectable by LOFAR. However, for 20 MeV LDM the 
situation is far more optimistic with significant values of 
5Tb — 5Tb,o at all times z < 800, reaching 35mK at z ~ 50, 
and approximately equal to 10 mK at z = 12, easily within 
the estimated sensitivity of LOFAR. 

An interesting extension of our analysis would be to 
consider other d ark matter particles, like exciting dark 
matter (XDM) jFinkbeiner fc Weinerl l2007t ). XDM can 
annihilate to produce two intermediate scalars <j> that 
can subsequently decay to standard model particles. If 
2m e < m$ < 2m M , the <f) will mainly decay to an e + e~ pair, 
with an energ y spectrum extending up to the mas s of the 
XDM particle (|Cholis. Goodenough fe Weinerl \200j ). 
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APPENDIX 

In this section we conveniently list the relevant parame- 
ters associated with each of the clumping factors utilised 
throughout this study, for clumping factors calculated us- 
ing structures possessing Moore (table 1) and NFW (table 2) 
profiles for our four neutralino dark matter models, and for 
Burkert profiles using our two LDM candidates (table 3). 
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"min. "cut F suh ,,F sa (%) N c Model 1 Model 2 Model 3 Model 4 



Ml 


-12 


-12 


3 


3 


N 




N 


N 


N 


M2 


-12 


-12 


3 


1.5 


N 




N 


N 


N 


M3 


-12 


-12 


0.3 


3 


N 




N 


N 


N 


M4 


-12 


-12 


0.3 


1.5 


N 




N 


N 


Y, (-0.12) 


M5 


-12 


6 


3 


3 


N 




N 


N 


N 


M6 


-12 


6 


3 


1.5 


N 




N 


N 


Y, (-0.04) 


M7 


-12 


6 


0.3 


3 


N 




N 


Y, (0.02) 


Y, (-0.04) 


M8 


-12 


6 


0.3 


1.5 


N 


Y. 


(-0.23) 


Y, (-0.08) 


Y, (-0.04) 


M9 


- 1 


- 1 


3 


3 


N 




N 


N 


N 


M10 


-4 


-4 


3 


1.5 


N 




N 


N 


N 


Mil 


-4 


-4 


0.3 


3 


N 




N 


N 


Y, (-0.07) 


M12 


-1 


-1 


0.3 


1.5 


N 




N 


N 


Y, (-0.06) 


M13 


-1 


6 


3 


3 


N 




N 


Y, (0.15) 


N 


M14 


-4 


6 


3 


1.5 


N 




N 


Y, (-0.21) 


Y, (-0.03) 


M15 


-4 


6 


0.3 


3 


N 




N 


Y, (-0.21) 


Y, (-0.03) 


M16 


-1 


6 


0.3 


3 


N 


Y. 


(-0.23) 


Y, (-0.09) 


Y, (-0.03) 


M17 


6 


6 


3 


3 


N 


Y. 


(-0.05) 


Y, (-0.09) 


Y, (0.00) 


M18 


6 


6 


3 


1.5 


Y, (-0.27) 


Y. 


(-0.05) 


Y, (-0.09) 


Y, (0.00) 


M19 


6 


6 


0.3 


3 


Y, (-0.27) 


Y, 


(-0.05) 


Y, (-0.09) 


Y, (0.00) 


M20 


6 


6 


0.3 


1.5 


Y, (-0.27) 


Y. 


(-0.05) 


Y, (-0.09) 


Y, (0.00) 



Table 1. Parameters for the clumping factors used to generate the results presented for structures possessing Moore density pro- 
files. For each clumping factor, the following parameters are displayed in each column (from left to right): column (1) - indicates the 
clumping factor designation; column (2) - n m j n =logio(M m j n /Mq), where M m j n is the minimum halo mass considered; column (3) 
- icut. =logio(Af cu t./MQ), where M cu t. is the mass below which the concentration-mass relation for halos is truncated; column (4) - 
percentage of host halo (subhalo) mass per mass decade in subhalos (sub-subhalos), column (5) - ratio of concentrations for a subhalo 
and halo of the same mass located at the same redshift; column (6) - indicates whether the relevant clumping factor satisfies the criterion 
discussed in §[7]regarding constraints on the diffuse background ('N' for No, 'Y' for Yes) for neutralino dark matter described by model 1; 
Column (7) (also under the heading 'model 1') - indicates the value of the differential brightness temperature difference relative to the 
standard "no dark matter" scenario, <5T(, — ST^ q (mK), evaluated at redshift z = 12 (provided that column (5) is 'Y') for neutralino dark 
matter described by model 1; Columns (8)&(9) are the same as for columns (6)&(7) but for neutralino dark matter described by model 2; 
Columns (10)&(11) are the same as for columns (6)&(7) but for neutralino dark matter described by model 3; Columns (12)&(13) are 
the same as for columns (6)&(7) but for neutralino dark matter described by model 4. 



n mill , n cut F Buh ,,F as (%) N c Model 1 Model 2 Model 3 Model 4 



Nl 


-12 


-12 


3 


3 


N 




N 




N 




N 


N2 


-12 


-12 


3 


1.5 


N 


Y, 


(-0.21) 


V, 


(-0.24) 


Y. 


(-0.02) 


N3 


-12 


-12 


0.3 


3 


N 


Y, 


(-0.18) 


V, 


(-0.23) 


Y. 


(-0.02) 


N4 


-12 


-12 


0.3 


1.5 


Y, (0.31) 


V. 


(-0.13) 


Y, 


(-0.20) 


Y 


-0.01) 


N5 


-12 


6 


3 


3 


N 




N 


V, 


(-0.24) 


Y. 


(-0.02) 


N6 


-12 


6 


3 


1.5 


Y, (-0.23) 


Y, 


(-0.05) 


V, 


(-0.08) 


Y 


(0.00) 


N7 


-12 


6 


0.3 


3 


Y, (-0.23) 


Y, 


(-0.05) 


Y 


(-0.08) 


Y 


(0.00) 


N8 


-12 


6 


0.3 


1.5 


Y, (-0.22) 


Y 


-0.04) 


Y 


(-0.07) 


Y 


(0.00) 


N9 


-4 


-i 


3 


3 


N 




N 




N 


Y 


-0.02) 


N10 


-4 


-4 


3 


1.5 


Y, (-0.15) 


Y, 


(-0.09) 


v. 


(-0.14) 


Y. 


(-0.01) 


Nil 


-4 


-4 


0.3 


3 


Y, (-0.17) 


Y, 


(-0.08) 


v. 


(-0.13) 


Y 


-0.01) 


N12 


-4 


-4 


0.3 


1.5 


Y, (-0.20) 


V. 


(-0.07) 


Y, 


(-0.11) 


Y 


-0.01) 


N13 


-1 


6 


3 


3 


N 


Y, 


(-0.07) 


Y 


(-0.11) 


Y. 


(-0.11) 


N14 


-4 


6 


3 


1.5 


Y, (-0.21) 


Y, 


(-0.04) 


Y 


(-0.06) 


Y 


(0.00) 


N15 


-4 


6 


0.3 


3 


Y, (-0.21) 


Y, 


(-0.04) 


Y 


(-0.06) 


Y 


(0.00) 


N16 


-1 


6 


0.3 


3 


Y, (-0.20) 


Y 


-0.03) 


Y 


(-0.06) 


Y 


(0.00) 


N17 


6 


6 


3 


3 


Y, (-0.04) 


Y 


(0.00) 


Y 


(-0.01) 


Y 


(0.00) 


N18 


6 


6 


3 


1.5 


Y (-0.04) 


Y 


(0.00) 


Y 


(-0.01) 


Y 


(0.00) 


N19 


6 


6 


0.3 


3 


Y (-0.04) 


Y 


(0.00) 


Y 


(-0.01) 


Y 


(0.00) 


N20 


6 


6 


0.3 


1.5 


Y (-0.04) 


Y 


(0.00) 


Y 


(-0.01) 


Y 


(0.00) 



Table 2. Same as for table 1 but for clumping factors associated with structures possessing NFW density profiles. 
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^min. 


Hcut 


F B ub. i Fss (%) 


N c 


wildm = 20 MeV 


wldm = 3 MeV 


Bl 


-0.80 


-0.80 


3 


3 


N 




132 


-0.80 


-0.80 


3 


1.5 


N 




133 


-0.80 


-0.80 


'3 


3 


N 




13 1 


-0.80 


-0.80 


0.3 


1.5 


N 




13 o 


-0.80 


C) 


3 


3 


N 




136 


-0.80 




3 


1.5 


N 




137 


-0.80 


g 


0.3 


3 


Y fl 1 03*1 

A . 111. KJrJ J 




138 


-0.80 


6 


0.3 


1.5 


Y, (10.70) 




139 


1.66 


1.66 


3 


3 




N 


BIO 


1.66 


1.66 


3 


1.5 


- 


N 


Bll 


1.66 


1.66 


0.3 


3 




N 


B12 


1.66 


1.66 


0.3, 


1.5 




Y, (-0.18) 


B13 


1.66 


6 


3 


3 




N 


B14 


1.66 


6 


3 


1.5 




N 


B15 


1.66 


6 


0.3 


3 




Y, (-0.15) 


B16 


1.66 


6 


0.3 


1.5 




Y, (-0.15) 


B17 


6 


6 


3 


3 


Y, (2.07) 


N 


B18 


6 


6 


3 


1.5 


Y, (2.04) 


Y, (-0.15) 


B19 


6 


6 


0.3 


3 


Y, (1.90) 


Y, (-0.04) 


B20 


6 


6 


0.3 


1.5 


Y, (1.90) 


Y, (-0.04) 



Table 3. Same as for table 1 but for clumping factors associated with structures possessing Burkert density profiles and using LDM of 
mass 20 MeV ( columns (6)&(7) ) and 3 MeV ( columns (8)&(9) ). 



